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We study nonlocal resistance in an H-shaped two-dimensional HgTe/CdTe quantum well consist 
of injector and detector, both of which can be tuned in the quantum spin Hall or metallic spin Hall 
regime. Because of strong spin-orbit interaction, there always exist spin Hall effect and the nonlocal 
resistance in HgTe/CdTe quantum well. We find that when both detector and injector are in the 
quantum spin Hall regime, the nonlocal resistance is quantized at 0.25^, which is robust against 
weak disorder scattering and small magnetic field. While beyond this regime, the nonlocal resistance 
decreases rapidly and will be strongly suppressed by disorder and magnetic field. In the presence of 
strong magnetic field, the quantum spin Hall regime will be switched into the quantum Hall regime 
and the nonlocal resistance will disappear. The nonlocal signal and its various manifestation in 
different hybrid regimes originate from the special band structure of HgTe/CdTe quantum well, and 
can be considered as the fingerprint of the helical quantum spin Hall edge states in two-dimensional 
topological insulator. 

PACS numbers: 72.20.-i, 73.20.-r, 72.25.Dc, 73.63.-b 


I. INTRODUCTION 

A topological insulator is a special quantum matteili^. 
Due to its particular band structure that the bulk states 
have a gap and the surface states can exist in the bulk 
band gap, the topological insulator behaves as an insu¬ 
lator in its interior and behaves as a metal on the sur¬ 
face. This leads to quantum spin Hall (QSH) effect. Dif¬ 
ferent from magnetic field induced quantum Hall effect 
where the time-reversal symmetry is broken, the QSH ef¬ 
fect arises from strong spin-orbit interaction (SOI) and 
is protected by the time-reversal symmetry. In the QSH 
regime, the electron spins are locked to their momenta 
and the boundary states are helical, i.e., there are two 
time reversed counter propagating edge states occupied 
by electrons with opposite spin. Such a topological state 
can be signed by Z 2 indesP. In topological insulator, the 
boundary states are protected by the time-reversal sym¬ 
metry and the backscattering between boundary states is 
strongly suppressed. As a result, the helical edge states 
are robust against the non-magnetic disorder.*^ In the 
bulk energy band beyond the QSH phase, the topologi¬ 
cal insulator behaves as a metal, which is different from 
the conventional two-dimensional (2D) metal. In con¬ 
ventional metal, the electron wave functions are local¬ 
ized as long as there exists any weak disordeP. How¬ 
ever, in topological insulator, the metal state can exhibit 
quantum conductance in the moderate disorder, which is 
called topological Anderson insulator phenomenePt^l. 

Up to now, the topological states have been predicted 
theoretically in several materials, such as the HgTe/CdTe 
quantum welP^^, the InAs/GaSb quantu m welP ^^*^, 
the monolayer graphene with intrinsic SOP^EU^ 
the gated bilayer graphen^Sl that contains the one¬ 
dimensional chiral edge states. In experiment, the 
HgTe/CdTe quantum well and the InAs/GaSb quantum 
well with inverted band have been successfully discov¬ 


ered as 2D topological insulators with the QSH phas^. 
Since the discovery of the topological insulator, many 
works have been concentrated on the verification of its 
helical edge states. Such as, Konig and co-workerP ob¬ 
served a quantized longitudinal conductance at about 
2e^/h that is consistent with the number of edge states 
predicted theoretically. Then, the transport along the 
edge states was confirmecP^E^. However, it remained to 
be shown that whether the transport due to the helical 
edge states is spin polarized. For this purpose, Br-iine and 
co-worker^i3 designed another experiment, in which the 
QSH effect and metallic spin Hall (MSH) effect are com¬ 
bined in a single HgTe/CdTe quantum well device using 
split gate technique. Through observation of nonlocal 
resistance in a H-shaped device, the spin polarization of 
edge state is then determined. In the work, in order to 
estimate the nonlocal resistance, the semiclassical sim¬ 
ulation are performecP^, but it breaks down when the 
chemical potential is close to the insulating gap. It is 
nowhere near enough for the QSH system. On the other 
hand, the nonlocal transport originates from the Hall ef¬ 
fect, including the quantum Hall effect and the spin Hall 
effect. So, we must carefully examine the nonlocal ef¬ 
fect to illustrate the role of the QSH state. Furthermore, 
as shown in the above experiments, the nonlocal signal 
deviated from the standard pattern predicted theoreti¬ 
cally because of the various impurity and dephasing ef¬ 
fect. Therefore, the detailed mechanism of the nonlocal 
transport in HgTe/CdTe quantum well, especially for the 
nonlocal transport in the QSH regime, is not very clear. 

In this paper, based on a four-band tight-binding 
model and aided by Keldysh nonequilibrium Green’s 
function, we study the nonlocal transport in a hybrid 
HgTe/CdTe quantum well, especially in the QSH regime. 
Following the recent experiment by Brzine et. al.“ we 
consider an H-shaped device based on the HgTe/CdTe 
quantum well as shown in Fig. 1(a). The on-site energies 
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FIG. 1: (Color online) (a) Sketch of the H-shaped hybrid 
four-terminal system. Through spilt gate technique, the on¬ 
site energies El and E2 in the injector (green region) and 
detector (yellow region) can be tuned separately. In the sys¬ 
tem, the scattering region (shadowed region) is connected to 
four leads, (b) Band structure of HgTe/CdTe: the conduction 
band (n-type carriers corresponding to n-MSH effect) and the 
valence band (p-type carriers corresponding to p-MSH effect) 
are separated by bulk band gap which holds the helical edge 
states and consequently the QSH effect. 


quantum well. The details will be expatiated in the sec¬ 
tion III. 

The rest of the paper is organized as follows. In Sec. II, 
based on the four-band tight-binding representation, the 
model Hamiltonian of system including central scattering 
region and attached leads is introduced. The formalisms 
for calculating the Green’s functions and nonlocal resis¬ 
tance are then derived. Sec. Ill gives numerical results 
together with detailed discussions. Finally, a brief sum¬ 
mary is presented in Sec. IV. 


II. MODEL AND HAMILTONIAN 

The whole system we consider is composed of the 
HgTe/CdTe quantum well with inverted band. In gen¬ 
eral, the structure of quantum well is asymmetry, which 
leads to the external Rashba SOI. The Hamiltonian of 
the system can be written as H{k) = Ho{k) + Hn^k), 
where Hji comes from the Rashba SOI, and 

H„(k) = ( ) . ( 1 ) 

From the time-reversal symmetry, we can get H^^(k) = 
H*{—'k), where 


El and E2 in the bottom (green) and top (yellow) re¬ 
gions can be tuned separately by the split gates above 
the two regions. As a result, the two regions can then be 
in the QSH or MSH regime [see Fig.1(b)]. When chang¬ 
ing El and i?2, the quantum well will be in different 
regime, i.e., different hybrid structure, including QSH- 
QSH regime, QSH-MSH regime, and MSH-MSH regime. 
When injecting current from the bottom two terminals, 
spin accumulation or spin polarized potential is gener¬ 
ated in the bridge between the bottom and top terminals 
due to the spin Hall effect. Then, because of the in¬ 
verse SHE, the charge voltage will be detected in the top 
two terminals and leads to the nonlocal resistance in the 
HgTe/CdTe quantum well. It is obvious that there ex¬ 
ists nonlocal resistance in all the hybrid devices because 
of the spin Hall effect. However, their manifestation is 
very different. For example, the nonlocal resistance in 
QSH-QSH device is quantized with the value of that 
is robust against moderate disorder and weak magnetic 
field, while the nonlocal resistance in MSH-QSH or MSH- 
MSH device is oscillating and fragile. In the presence of 
strong magnetic field, the quantized nonlocal resistance 
will shrink and disappear finally. All these features of the 
nonlocal resistance: the robust quantized nonlocal resis¬ 
tance in the QSH regime, the fragile nonlocal resistance 
in the MSH regime, and the vanishing nonlocal resistance 
in the presence of strong magnetic held, are the partic¬ 
ular characters of HgTe/CdTe quantum well, which can 
be derived from the special band structure of HgTe/CdTe 


H^{k) = hvFikxCTx -kyay) + {m + Ck ■'k)az 

-)- Dk ■ k(TQ, (2) 

where crx,y,z are Pauli matrices presenting the pseudo 
spin formed by s and p orbitals, (Tq is the 2x2 unitary 
matrix in the pseudo spin space. Ui? is Fermi velocity. 
C, D, and m are the system parameters, which can be 
experimentally controlled. While the Hamiltonian due 
to Rashba SOI Hu is expressed as 
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where tn is the Rashba SOI strength. 

Equations ([^ and (§ are the low-energy effective 
Hamiltonians of the HgTe/CdTe quantum welPE^ from 
the k -p perturbation. For an ideal crystal lattice, i.e., 
the inhnite periodic system, the momentum fc is a good 
quantum number. In this case, using Hamiltonians Q 
and (§ is convenient. However, here we consider charge 
transport in the H-shaped system, Hamiltonian can’t be 
expressed in momentum space, it should be expressed 
in real space. To do this, we substitute k with —IV 
and use the finite-difference approximation, then the to¬ 
tal Hamiltonian El is transformed into the tight-binding 
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Hamiltonian in square lattice. It is given 
H = Y^dlHndi 

i 

+ ^ 46*'^“'*+“- di+a, + h.C. 
i 

,i+aj,di+aj, + h.C., (4) 

i 

where di = [ds,i,t) dp,i.;]^ with ‘T’ denoting 

transpose, and ds(p) 4 ,^(^) and j are the annihila¬ 
tion and creation operators for s or p orbital at site i with 
spin up or spin down, respectively, i = (ixiiy) is the in¬ 
dex of the discrete site of the system in the square lattice, 
and = (a, 0 ) and a.y = ( 0 , a) are the unit vectors of the 
square lattice with a the lattice constant. In zero mag¬ 
netic field, the Hamiltonian @ possesses time-reversal 
invariant. In the presence of a uniform perpendicular 
magnetic field B = [0,0,5], the time-reversal symme¬ 
try is broken. In coulomb gauge, considering the seimi- 
infinite leads being along the cc-direction, the vector po¬ 
tential is set as A = [—By, 0, 0] which is y dependent 
but periodic in the ^-direction. In this case, a phase 
0 ij is generated in the hopping term idj ■ The phase 
d’ij = /^ A • d\/(j)Q with flux quanta (j)Q = h/e. In Eq.(|^, 
Hii and idi,i+a,„(a„) are all 4 x 4 block matrix that are 
expressed as 


Ha = (ci - ^)(so (8) (To) + {m- 4-^)(so 8 cFz) 

rr D . , C , 

Hi.i+e,,, = ^(so 8 (To) + ^(So 8 (Tz) 

f „ , , .tR f „(To + (Tz, 

-*^(5.8(t.)+*-(.,8^—) 

zr D , , C , 

= ^(so 8 (To) + ^(so 8 (Tj,) 

, -^^F , „ X .tR, „(To + (Tz, , x 


where Sx,y,z are the pauli matrices denoting the real spin 
and So is the 2 x 2 unitary matrix extented in real spin 
space. €i = E1{E2)+Wi with El and E2 being the on-site 
energies in detector (yellow) and injector (green) regions 
in Fig. 1(a). Wi comes from the disorder effect that is a 
random on-site potential which is uniformly distributed 
in the region [—w/ 2 ,w/ 2 ]. 

Based on above Hamiltonian, the charge current flow¬ 
ing to the p-th lead can be calculated from the zero tem¬ 
perature Landauer-Buttiker formulcP^ 

Jp='^J 2 ^P,{Vy-V,) ( 6 ) 

where p,q = 1 , 2 ,3,4 are the index of the four leads. 
Tpq is the transmission coefficient from terminal q to 
terminal p. In the following, we derive Tpq in the sys¬ 
tem without and with Rashba SOI. In the absence of 


Rashba SOI, the Hamiltonian with spin up and spin 
down are decoupled, and the self energy and Green’s 
functions for the system with spin up and spin down 
can be calculated separately. Then transmission Tpq = 
Tpq.t + Tpq,(y{E) = Trr^o[^p,aGpq,^Tq^^Gqp,^], 

where the ”Tiv,o” denotes the trace over real space and 
orbital space (s and p orbitals), the line-width function 
Tg.cr = ~ Gpq ^ is the Green’s function 

matrix whose rows and columns mark the lattices that 
are nearest to p and q lead respectively. The Green’s 
function G;(5) = G“;t(F) = (El - H, -Eq^.a)-\ 
where Hu is Hamiltonian matrix with spin cr in the cen¬ 
tral region and I is the unit matrix with the same di¬ 
mension as that of Hu, and is the retarded self¬ 
energy function contributed by the electrons with spin 
(T in lead q. The self energy function can be obtained 
from YTp u = Hcp.ugp^uHpc.a, where Hcp,a is the coupling 
from central region to lead p and Pp u is the surface re¬ 
tarded Green’s function of semi-infinite lea d p wh ich can 
be calculated using transfer-matrix methocP^lSU When 
considering Rashba SOI, electrons spins (spin up and 
spin down) are coupled with each other and the total 
transmission Tpq{E) = Trr^o.sppGp^jTqG^p], where the 
” Tr,.,o,s” is the trace over real space, orbital space and 
spin space. 

In the following, we calculate the nonlocal response, 
i.e., the voltage response (detector) in top region on the 
current (injector) in bottom region, which can be denoted 
by nonlocal resistance i? 23 ,i 4 - We also calculator i?i 4 , 23 j 
i.e., the voltage response in bottom region on the current 
in top region. i? 23 ,i 4 and i?i 4,23 have the nearly same 
properties, so, we focus on i? 23 ,i 4 in the following. We 
apply a bias V across the injector terminals 1 and 4 to 
inject current. For the detector terminals 2 and 3, the 
currents are set to zero. Then use the boundary condi¬ 
tions Vi = V, V 4 . = 0 , J 2 = J 3 = 0 , we can calculate the 
current Ji = — J 4 and the voltages V 2 and V 3 in the volt¬ 
age probes from the Landauer-Buttiker formula. Finally, 

the nonlocal resistance is given by i?23,14 = {V 2 —Vi)/Ji. 


III. NUMERICAL RESULTS AND ANALYSIS 

In the numerical calculation, the parameters of 
the HgTe/GdTe quantum well are set as hvp = 
364.SmeUnm, G = — 686 meVnm^, D = —512TOeynTO^, 
and the effective mass is taken as m = — lOmeU, 
which corresponds to the realistic quantum well with 
thickness d = 7nm^ It exceeds the critical thickness 
dc = 6.3nm and induces the inverted band, which leads 
to the topological phase. Gomparing to the inverted 
InAs/GaSb/AlSb quantum well, the Rashba SOI in 
HgTe/GdTe quantum well is very small and can be usu¬ 
ally neglected in the numerical calculation. However, in 
order to quantitatively estimate the effect of Rashba SOI, 
we set a nonzero Rashba SOI. The strength of Rashba 
SOI is set to a = 50meVnm. It is a very large value, 
because the Rashba SOI in inverted InAs/GaSb/AlSb 
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FIG. 2: i? 23 ,i 4 vs on-site energy El in the injector region 

with Rashba SOI a = 0 (black solid lines) and SOmeVnm 
(red dotted lines). Panels (a), (b), (c), and (d) correspond to 
E2 = — 12meF, —5meV —2>QmeV and 2bmeV , respectively, 
where E2 is the on-site energy in the detector (top region). 
The gray region in panel (a,b,c,d) is the range of bulk energy 
gap of injector (the sample bottom). 


quantum well is only a = TlmeVnmp^ while Rashba 
SOI in HgTe/CdTe quantum well is much smaller than 
that in InAs/GaSb/AlSb quantum well. Furthermore, 
we set lattice constant a = 5nm, that is small enough to 
get a reasonable band structure. The scattering region 
is shadowed in Fig. 1(a). The width and length of scat¬ 
tering region is set to W = Na = 750nm with N = 150 
and L = Ma = 750nm with M = 150. As shown in 
Fig. 1(a), the width of the bridge that connects the top 
and bottom terminals, and the distance between top and 
bottom terminals, are all 250nm. In the numerical cal¬ 
culation, we fix the Fermi energy at Ep = 0 and change 
the on-site energies El and E2. In fact, we can also fix 
El = E2 = 0 and tuning the Fermi energy in the two 
regions. They are equal in the calculation. 

Based on this device, we first study the nonlocal re¬ 
sistance in the clean hybrid system without the external 
magnetic field. With the change of the on-site energies 
El and E2 in the injector and the detector, there will be 
three different hybrid regimes: QSH-QSH, QSH-MSH, 
MSH-MSH. The characters of the nonlocal response in 
these hybrid regimes are depicted in Fig.2. In panel 
(a),(b),(c) and (d), we plot the i? 23 ,i 4 vs El in the bot¬ 
tom region (injector) with and without Rashba SOI for 
E2 = —12meV, —5meV, —30meV and 25meV, respec¬ 
tively. Here, in the calculation, we set the Fermi energy 
Ep = 0. And E2 = —30meV, —12meV, —5meV, and 
25meV correspond to the the detector (top region) is in 
the n-MSH, near QSH, QSH, and p-MSH regimes, re¬ 
spectively. Similarly, when El changes from —40meV 
to 30meR, the injector (bottom region) develops from n- 
MSH to p-MSH regime via QSH regime. Now, we analyse 



FIG. 3: (Golor online) R 23 ,i 4 vs on-site energy El in the 
injector region for different disorder strength w = 20meV, 
50meV, lOOmeF, and 200meV. Panels (a), (b), (c), and (d) 
correspond to E2 — —12meV, —5meV, —30meV and 25meV, 
respectively. 


the nonlocal properties in different hybrid structure. 

From Fig.2, we can find some common characters: 
(1) The Rashba SOI hardly affects the nonlocal effects 
of hybrid HgTe/CdTe quantum well in variant regimes 
(QSH-QSH, QSH-MSH, MSH-MSH regimes, and so on). 
For the system with a = 0 (black solid lines) and 
a = 50meVnm (red dotted lines), the nonlocal resistance 
are almost the same, there are only slightly quantitative 
difference between them. So the Rashba SOI is unim¬ 
portant on studying the nonlocal effect of topological in¬ 
sulator. (2) No matter what the values of E2^ R 23 ,i 4 is 
always biggest when the injector is in the QSH regime, 
i.e., El G [—10meV,10meV]. This means the nonlocal 
effect is most remarkable in the QSH regime, that is well 
consistent with the experiment in Ref. |17j. (3) i? 23 ,i 4 
in the n-MSH and p-MSH regimes is small and shows 
oscillating behavior, it can even be negative for some 
special El. Furthermore, the value of i? 23 ,i 4 in the p- 
MSH regime is bigger than that in the n-MSH regime, 
as shown in the experiment.^ (4) At El « —7.5meV, 
because of the finite size effect, i? 24 .i 3 abruptly peaks or 
dips. Now we analysis Fig.2 in detail. We focus on several 
main regime: the QSH-QSH regime [the central region of 
Fig.2(b)], the QSH-MSH regime [the left and right re¬ 
gions of Fig.2(b)], and the MSH-MSH regime [the left or 
right panel in Figs.2(c) and 2(d)]. When the injector and 
detector regions are all in the QSH regime (QSH-QSH hy¬ 
brid regime), because of the counter propagating helical 
edge states, i? 23 .i 4 is biggest and is quantized at the value 
of regardless of the Rashba SOI [see the central re¬ 
gion in Fig.2(b)]. It is interpreted as follows. In the pres¬ 
ence of the helical edge states, the transmission coeffi¬ 
cients are integer T 12/21 = 2 ) 23/32 = 234/43 = ^ 41/14 = 1 , 
we can then conclude V 2 = -I- |14, V 3 = ^Vi + |V 4 , 
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FIG. 4: Nonlocal resistance i? 23 ,i 4 vs on-site energy El 
in the injector region. The on-site energy in the detector 
region is set as E2 = —bmeV. Panels (a), (b), (c), and (d) 
correspond to the magnetic field B = 0, O.OST, 0.16T and 
0.32T, respectively. 


and Ji = —J 4 = (4/3)(Fi — V 4 ) from the Landauer- 
Buttiker formula. Then, i? 23 ,i 4 = = 0.25/i/e^. 

Next we consider that the energy level of the detector is 
E2 = —12meV. In this case, although the Fermi energy 
of the detector is in the conduction band, it is very near 
the band edge, lOmeV, in which the bulk states coexist 
with edge states that leads to some exotic behaviors. We 
will call it the near QSH regime. When the detector is 
in the near QSH regime and the injector is in the QSH 
regime, the hybrid system will be in the near-QSH-QSH 
regime. In this regime, i? 23 .i 4 is decreased by the bulk 
states. However, the nonlocal effect is still strong [see 
Fig.2(a)] due to the helical edge states that survive in 
the bulk near band edge. Beyond the QSH-QSH regime, 
the MSH effects paly a role on the nonlocal transport 
and i? 23 ,i 4 decreases because the spin Hall effect in MSH 
phase is weaker than that in QSH phase. In the QSH- 
nMSH and QSH-pMSH regimes [see left and right re¬ 
gions of Fig.2(b)], because the eigenstates of the MSH 
system are extended in the transverse directiorP^, i? 23 ,i 4 
decreases rapidly and oscillates around zero. The oscil¬ 
lating frequency is coincident to the sub-band distribu¬ 
tion in the conduction and valence bands. Finally, in the 
MSH-MSH regime [Fig.2(c,d)], the nonlocal resistance is 
induced by MSH completely, so i? 23 .i 4 is oscillating and 
very small. In this case, the nonlocal resistance i? 23 ,i 4 
can be negative when El is of some special values. 


Since the Rashba SOI is not important in the nonlocal 
effect, only the system with a = 0 is studied in the fol¬ 
lowing. In Fig.3, we study the disorder effect of the non¬ 
local resistance R 23.14 in zero magnetic field. We choose 
an H-shaped scattering regions [the shadowed region in 
Fig. I (a)], in which the on-site potential Wi is randomly 
distributed in the region of [—f,f], with w the disor¬ 
der strength. When the detector (related to E2) is in 
the QSH regime [Fig.3(b)], it is possible for the system 
to be in the QSH-QSH, MSH-QSH and near-QSH-QSH 
regimes. In the presence of weak disorder {w = 20meV), 
R 23 ,i 4 is quantized in the QSH-QSH regime, and the os¬ 
cillating of i? 23 ,i 4 in the MSH-QSH regime is strongly 
depressed [see the solid lines in Fig.2(b) and Fig.3(b)]. 
In the moderate disorder (w = fiOmeH), i? 23 ,i 4 in the 
QSH-QSH regime is still kept but i? 23 .i 4 in QSH-MSH is 
increased [see red dashed line in Fig.3(b)]. Besides, due 
to the topological Anderson insulator phenomenorP^, 
R 23 M in the near-QSH-QSH regime, where the detec¬ 
tor is near QSH regime and the injector is in the QSH 
regime, is quantized by the moderate disorder [see the 
red dashed line in Fig.3(a)]. It means detector that is 
near QSH regime is now driven into the QSH regime by 
the moderate disorder. Then the disorder effect in QSH- 
QSH and near-QSH-QSH regime are almost same [see 
blue dotted lines in Fig.3(a) and Fig.3(b)]. When both 
El and E2 are far away from the energy gap, the system 
is in the MSH-MSH regime and i? 23 ,i 4 is increased by the 
weak disorder [see Fig.3(d)]. When the disorder is very 
strong, random scattering dominates the nonlocal trans¬ 
port and i? 23 .i 4 in all the regimes are nearly the same 
[see the green dash dotted lines Fig.3(a,b,c)]. 

Now, we study the nonlocal transport under the exter¬ 
nal magnetic field. We first consider the weak magnetic 
field. In the HgTe/CdTe quantum well, we are interested 
in the nonlocal properties of the helical edge states. So, 
in the following, we consider only the hybrid system with 
E2 = —bmeV, i.e., the detector is in the QSH regime. 
Fig.4 depicts i? 23 ,i 4 vs El for different weak magnetic 
field B. The corresponding band structures are plotted in 
Fig.5. From Fig.4, we can find with increasing magnetic 
field B, the quantized R 23.14 in the QSH-QSH regime 
hardly changes and i? 23 ,i 4 in the QSH-nMSH regime 
is suppressed gradually. This can be explained as fol¬ 
lows. As we know, the quantized i? 23 .i 4 in the QSH-QSH 
regime arises from the helical edge states, while i? 23 ,i 4 
in the QSH-nMSH regime is dominated by the extended 
eigenstateP^ in the MSH system. Considering the trans¬ 
port process, the helical edge states propagate only along 
the geometric edge, while the extended MSH states can 
propagate in the whole scattering region. Thus, when the 
external magnetic field is added (in all regions including 
the four leads and the scattering region), it can drasti¬ 
cally affect i? 23 ,i 4 in the QSH-MSH regime but hardly 
changes the quantized i? 23 ,i 4 - All in all, in the presence 
of weak magnetic field, it is the helical edge states and the 
extended MSH states that induce the unchanged quan¬ 
tized i? 23 ,i 4 in the QSH-QSH regime and drastically sup- 
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FIG. 5: (Color online) Band structure of HgTe/CdTe quantum well in the presence of weak magnetic field. The bottom and 
top panels are for spin up and spin down, respectively. Two gray lines denote the edge of band gap in the absence of magnetic 
field. The on-site energy is set to zero. 
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FIG. 6: (Color online) Nonlocal resistance i? 23 ,i 4 vs on-site 
energy El in the injector region for different magnetic field. 
The on-site energy in detector region is set as E2 = —5meV. 



press i? 23 ,i 4 in the QSH-nMSH regime, as in the weakly 
disordered system. Besides, in Figs.4(b), 4(c), and 4(d), 
we can see that in some energy regions, e.g., El near 
— lOmeV as in Fig.4(b), El G [—I5meV) —lOmeF] as in 
Fig.4(c) and El € [—25meF, —lOmeF] as in Fig.4(d), ex¬ 
cept some abrupt peaks, i? 23 ,i 4 is zero. This can be inter¬ 
preted as follows. When the magnetic field increases, the 
flat Landau level forms gradually in the conduction band. 
The topological edge state near the conduction band gap 
with spin down is first broken by the magnetic field, be¬ 


cause its chirality, i.e., the rotation direction along the 
scatter edge, which can be characterized by Chern num- 
berdUlSr 

is opposite to the edge state induced by mag¬ 
netic field. The opposite chirality also induces the slight 
dip in Landau levels [see the top panels in Fig.5], which 
leads to the abrupt peaks in the conduction bancP^ as 
shown in Fig.4(c) and (d). On the other hand, for the 
carriers with spin up, the helical edge state near topolog¬ 
ical band gap is kept [see bottom panels in Fig.5] because 
its chirality is same to the Landau edge state. Combin¬ 
ing the Landau gap in the system with spin down and 
the topological edge state in the system with spin up, 
the system is equal to a quantum Hall system, in which 
the edge state is unidirectional, since the edge state con¬ 
tributed by spin down is broken by the magnetic field. As 
a result, the chemical potential V 2 and V 3 of the detector 
terminals are determined by one of the adjacent source 
terminals. It means the chemical potential V 2 = Vi and 
V 3 = V 2 , which leads to zero nonlocal resistance i? 23 ,i 4 - 
Finally, there are also abrupt dips in quantized i? 23 ,i 4 ) as 
interpreted in the zero magnetic field case, which comes 
from the finite size effect. 

Next, we consider the nonlocal effect in the strong mag¬ 
netic field. In this case, the edge states for spin up and 
spin down all develop into flat Landau levels. Fig .6 plots 
the nonlocal resistance i? 23 ,i 4 vs on-site energy El in 
strong magnetic field B. Three characters are found: 
(1) Although i? 23 ,i 4 is quantized in the bulk gap of [- 
10meV,10meV], the quantized range gradually shrinks 
with increasing of B. (2) Except several abrupt peaks. 
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FIG. 7: (Color online) Band structure of HgTe/CdTe quantum well in the presence of strong magnetic field. The bottom and 
top panels are for spin up and spin down, respectively. The on-site energy is set to zero. The red lines denote the position of 
Fermi level relative to the on-site energy. 


7?23.i 4 is nearly zero in the whole conduction band. (3) 
7?23 .i 4 in the valence band is zero in some region, this re¬ 
gion expands with increasing of B. To understand these 
behaviors, we need to analysis the band influence induced 
by strong magnetic field. 

As we know in the HgTe/CdTe topological insulator, 
the chirality of the edge state induced by the inverted 
band is identical for electron and hole but opposite for 
carriers with spin up and spin down. On the other hand, 
in the strong magnetic field, the chirality of Landau edge 
state induced by the magnetic field is same for carriers 
with spin up and spin down but opposite for electron 
and hole. For all four classes of carrier, i.e., electrons 
and holes with spin up and spin down (signed as e f, e 
/if, h 1), when the two type edge states,i.e., edge state 
induced by inverted band and magnetic field (signed by 
‘E’ and ‘B’) have same chirality, they can coexist and 
boost. Otherwise, they are annihilated with each other. 
In our system, concerning the chirality of two types of 
edge states, Ch-\^E = Ce^,E = Ce^,B = Ce^,B = 1, 
Ch^^B = ChiM = Chi,E = Cei,E = -1, the edge states of 
e t and h I are strengthened, and those of h f and e I are 
destroyed in the strong magnetic field, as shown in Fig.7. 
However, since the edge states induced by the inverted 
band is more robust in the energy close to the conduc¬ 
tion band, the destroy is more gentle (see top panels in 
Fig.7). So, in the strong magnetic field, the topologi¬ 
cal edge states induced by inverted band are gradually 
destroyed in the region of if S [—lOmeF, 0] and nearly 
kept in the region of if S [0, lOmeV] as shown in Fig.7. 


And in the region of iff € [0, lOmeF] in Fig.6 (corre¬ 
sponding to if S [—lOmeF, 0] in Fig.7), one of two heli¬ 
cal edge states gradually disappears and the system can 
be regarded as a quantum Hall system with i? 23 ,i 4 being 
zero. Then we demonstrate the first character in the last 
paragraph. Furthermore, in the presence of strong mag¬ 
netic field, the Landau levels are completely formed in 
the conduction band. In the deep conduction band, the 
topological edge state induced by inverted band does not 
work and the unidirectional Landau edge states dominate 
the transport procession, leading to the zero ii 23 ,i 4 in the 
conduction band of ifl < — lOmeF. Besides zero i? 23 ,i 4 , 
there are also some abrupt peaks in the conduction band, 
it is because of the slightly dip in Landau level. This is 
for the second character. Finally, in the valence band, 
although the edge states of spin up are destroyed by the 
magnetic field, Landau levels have not been completely 
formed and there still exist extend states between Lan¬ 
dau levels, which leads to the nonzero i? 23 ,i 4 in valence 
band. Except these extent states, Landau levels are en¬ 
tirely gapped and i? 23 ,i 4 is then zero. This interprets the 
third character in Fig.6. 

Besides the three characters depicted in Fig.6, we can 
also expect when B is large enough, the edge states of h f 
and e I are also completely destroyed, and the topolog¬ 
ical gap [—10meV,10meV] in these system will become 
the real g^in which both the bulk and edge state are 
all absent.^ Then, the quantized region of i? 23 ,i 4 will 
disappear entirely. Furthermore, when B is very large, 
in the conduction and valence bands, the Landau edge 















































































IV. CONCLUSION 



disorder (meV) 

FIG. 8: (Color online) Nonlocal resistance i? 23 ,i 4 vs disorder 
strength w for different magnetic field. The on-site energies 
of the injector and detector regions are set as El = E2 = 
—5meV. 

state will dominate the transport and the system come 
into the quantum Hall regime, in which i? 23 ,i 4 is strictly 
zero in any hybrid regime. So, the nonlocal resistance 
in the H-shaped hybrid HgTe/CdTe quantum well will 
completely disappear in very large magnetic field. 

In the following, we study the disorder effect on the 
nonlocal resistance in the QSH-QSH regime in the pres¬ 
ence of strong magnetic field. We set the on-site energies 
in the injector and the detector as El = E2 = —5meV 
and Fermi energy Ep = 0, which equal to Ep = 5 and 
El = E2 = 0. With these parameters, helical QSH 
edge state can be maintained at very large magnetic field 
B = 6.4T [see Fig.7(dl) and (d2) in which Ep = 5 and 
E = 0], However, because the magnetic field destroys the 
edge state in the band edge, it weakens the robustness 
of the helical edge state, especially for the edge states 
located near the gap edge, as shown in Fig.8, in which 
.^ 23,14 vs disorder strength for different magnetic field 
is plotted. It can be seen when B = 6.4T, the quan¬ 
tized i? 23 ,i 4 is decreased at small disorder strength of 
w « 40meU, while for B = 0, 0.8T and 3.2T, quan¬ 
tized i? 23 ,i 4 can be maintained even at w = lOOmeU. 
This is because at B = GAT [Fig.7(dl) and (d2)], the 
helical band gap is violently shrunk into the region of 
[2meV,8meV] and Ep = 5meV is very close to the gap 
edge. For the magnetic held B < 3.2T, the band gap is 
wider of about [—5meV, IQmeV]. In this case, the Fermi 
energy is far away from the gap edge and the nonlocal 
resistance can be maintained in the strong disorder scat¬ 
tering. 


In summary, we have investigated the nonlocal trans¬ 
port in a H-shaped hybrid HgTe/CdTe quantum well. 
Three hybrid regime are considered, QSH-QSH regime, 
QSH-MSH regime and MSH-MSH regime. It is found in 
the QSH-QSH regime, the spin polarized edge states de¬ 
nominate the transport procession, the nonlocal effect is 
most strong, and the nonlocal resistance i? 23 ,i 4 is quan¬ 
tized in the value of 0.25^. While in the QSH-MSH 
device, due to the extended states in MSH effect, the 
nonlocal resistance is oscillating and much smaller than 
in the QSH-QSH device. In the MSH-MSH regime, the 
nonlocal effect nearly disappear, i? 23 ,i 4 can only reach 
the order of 0.01^. In the presence of the disorder, the 
quantized nonlocal resistance in the QSH-QSH device is 
robust because of the time reversal protected edge states. 
Near the QSH-QSH regime, the nonlocal resistance is 
enhanced and quantized due to the topological Ander¬ 
son insulator phenomenon. While in the QSH-MSH and 
MSH-MSH device, the oscillated nonlocal resistance is 
strongly restricted by disorder. Finally, the magnetic 
field effect is investigated. It is found the quantized non¬ 
local resistance in QSH-QSH regime can’t be affect by 
weak magnetic field, but the nonlocal resistance of QSH- 
MSH device is smoothed out by B. With the increasing 
of B, the Landau level forms, counter propagating edge 
states are replaced by chiral edge state, so the region of 
quantized nonlocal resistance decrease and disappear fi¬ 
nally. All these aforementioned nonlocal features can be 
derived from the special band structure of HgTe/CdTe 
quantum well. It is the unique property of HgTe/CdTe 
quantum well in topological insulator phase, it can be re¬ 
garded as the fingerprint of the helical QSH edge states 
in 2D topological insulator. 
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